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Abstract: We use holographic techniques to calculate the hrst thermal correction to 
the entanglement entropy of a cap-like region of a CFT dehned on a sphere, successfully 
reproducing the held theory result. Since this is an order-one correction to the entropy 
in the large-expansion, quantum corrections to the holographic entanglement entropy 
formula are essential. The bulk calculation is made tractable using the same technical 
machinery recently used to derive the linearized Einstein equations in the bulk from 
the hrst law of entanglement entropy in the CFT. 
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1 Introduction 


Herzog [1] has recently calculated the hrst thermal correcton to the entanglement en¬ 
tropy of a cap-like region B for a CFT dehned on i? x S^~^. The CFT is assumed to 
have a mass gap A and degeneracy g in the hrst excited state, in terms of which we 
have 

\T) = Sb{T) - Sb{ 0) = gAe-^^lM + o(e-^^) (1.1) 
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and the CFT is at 


1 Jo 

the opening angle of the cap is 200 , fid-i is the volume of 5''^“^, 
temperature T = We are working in units where the sphere has radius one. 

The result (1.1) is independent of N in the sense of large-iV CFTs, and so in the 
large-iV expansion should appear at 0(1). To calculate 0(1) corrections to the entropy 
holographically, we need to use the quantum generalization of the Ryu-Takayanagi [3] 
entropy formula due to Faulkner, Lewkowycz, and Maldacena [4]: 
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^Herzog argues in Ref. [1] that the expression on the right-hand side of (1.1) can shift depending 
on contributions of boundary terms, such as those coming from the conformal coupling of a free scalar 
to the background curvature. In Ref. [2], the claim is that a generic interacting CFT will not have 
such boundary terms. We will reproduce (1.1) as written, and leave to future work the incorporation 
of boundary terms for non-generic interacting theories. 
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where A is the area of the extremal surface anchored on the boundary of B and S'buik is 
the entanglement entropy of the bulk matter in the region between B and the extremal 
surface.^ We will also note that, in hnding the difference between the T = 0 and hnite 
T entropies, there are some fruitful cancelations. From the CFT point of view there 
are short-distance divergences in the entropy which cancel in the difference; in the 
bulk this comes from cancelations near the boundary at inhnity. But there are similar 
short-distance divergences in S'buik that also cancel, which we would otherwise have to 
regulate with counterterms. To summarize, we must compute 

+ SStu,u(T), (1.4) 

4GrAr 

and show that it is equal to where the vacuum-subtracted area, SA(T), and 

vacuum-subtracted bulk entropy, dSbuik(T), are both hnite, well-behaved quantities. 

At high temperatures, above the Hawking-Page phase transition, we normally say 
that the geometry dual to the thermal state is a large black hole. Alternatively, one can 
say that all of the typical pure states in the canonical ensemble at high temperature 
look like the same large black hole, and so thermal expectation values can be computed 
using that black hole background.^ Below the Hawking-Page phase transition, this is 
not the case. In the low-temperature limit, the thermal state can be expanded as 

p(T) = |0)(0| + (|A)(A| - |0)(0|) + ■ ■ ■ . (1.6) 

Here we are denoting by |0) the vacuum state (i.e., empty AdS), and |A) the bulk held 
state with energy A (for notational simplicity we are assuming no degeneracy), which 
for example could be the lowest-energy single-particle state of a free scalar held of mass 
= A(A — d). The bulk represented by this mixed state not a single geometry, 
but is just what the equation implies: an incoherent mixture of the bulk vacuum and 
the lowest-lying excited bulk state. ^ We need to compute the extremal surface area in 
this mixed state, as well as the bulk entanglement entropy, and subtract the respective 
vacuum values. We will see that both of these computations reduce to calculations in 
the state |A), which we will perform. 

^For general theories of matter coupled to gravity, we would have to use a Wald-like entropy formula 
in place of the area and the associated terms. While we restrict ourselves to theories where entropy is 
represented by area for simplicity, the techniques of Ref. [5] that we employ, based on the formalism 
of Iyer and Wald [6] , can be used in the more general case. 

^Up to subtleties involving the near-horizon region. 

^The “classical” 0(1) part of the geometry is of course the same in both cases, but for our purposes 
the “quantum” 0(1/N^) part is also needed. 
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Figure 1. (Left) The r = 0 slice of AdS. The region B is on the boundary, B is the bulk 
extremal surface, and S is between them in the bulk. (Right) A cross-section of AdS, where 
time runs vertically. S lies at r = 0 and the causal wedge associated to B is shown. 



2 Geometric Setup 

The region B on the sphere whose entropy we are computing is dehned by its opening 
angle 26q. In terms of the polar angle 6*, it is the region 6 < 6q. In the bulk, we will 
typically use global coordinates where the metric takes the form 

(—dr'^ + dp^ + sin^ pdfl? i) , (2.1) 

cos^'p ^ ^ 

where dQ‘^_^ is the metric of a unit (d — l)-sphere, and the angle 6 is one of the 
coordinates on this sphere: 

d^a-i = ddQ‘^_2. ( 2 . 2 ) 

The boundary is located at p = ti/2. The extremal surface B which shares its boundary 
with B and is used to calculate the entropy is given by the equation 

sin p cos 9 = cos 9q. (2.3) 

The region between B and B on the r = 0 surface will be denoted by S (see Fig. 1). 
When we compute the bulk entanglement entropy, we can think of it as computing the 
entanglement entropy of the state restricted to S. 

It so happens that, since we are considering such a simple surface in empty AdS, 
B is also the boundary of the “causal wedge” associated to B. The causal wedge in the 
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bulk is the set of all points whose past and future intersects the domain of dependence 
of B on the boundary. We can cover the causal wedge with coordinates t, C, and u 
dehned by 


coshCsinhu = tan p sin 6^, 

• 1 A • 1 sinr 

sinh C sinh t = -, 

cos p 

, , , n ■ ^ ^ ^ ■ n COST 

cosh C cosh u + cos oq sinh Q cosh t = sin b'o-, 

cos p 

cos 6q cosh C cosh u + sinh ^ cosh t = sin 9q tan p cos 0, 
in terms of which the metric is of AdS-Rindler form: 


(2.4) 

(2.5) 

(2.6) 
(2.7) 


ds^ = — sinh^ C,dt^ + d(‘^ + cosh^ C{du‘^ + sinh^ udVt^^_ 2 ). (2.8) 


We see that on the boundary, oo, this choice of coordinates induces an i? x 
geometry on the domain of dependence of B. The surface B is obtained by taking 
C —> 0 at hxed t. The Killing vector = dt generates time translations in this AdS- 
Rindler patch, and for calculations it is useful to have an expression for it in terms of 
the global coordinates: 


cos r sin p cos 6 — cos On ^ cos 6 cos p sin r ^ sin 6 sin r ^ 
r =-- -dr +-- dp - - - —de. 


sin On 


sin On 


sin p sin Oq 


(2.9) 


Notice that vanishes on B. We have normalized so that it has unit surface gravity, 
K = 1, where on the Rindler horizon. 


3 Holographic Entropy Calculation 

The Area Term 

First we will compute the area term. According to the Ryu-Takayanagi [3] prescription, 
we need to hnd the area of the extremal surface in the bulk which is anchored on 
the boundary of R. As is typically formulated, this term only makes sense when 
there is a well-dehned notion of a classical geometry in the bulk. For our purposes 
a slight generalization is necessary. We will interpret the Ryu-Takayanagi term in 
the entropy to be the expectation value of the operator A which encodes the area of 
the extremal surface. This operator is a functional of the metric operator, which for 
small perturbations we can write g + h with g the AdS vacuum metric. In the low- 
temperature limit, we use (1.5) to write the difference between the hnite-temperature 
area and the vacuum area as 

5A{T) = Tr(ip(T)) - (0| i |0) = ((A| A |A) - (0| A |0)) + ■ ■ ■ . (3.1) 
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We are only interested in evaluating this difference to leading order in G^v, and so we 
can expand (A| A |A) to first order in h. Up to operator ordering ambiguities (which 
do not appear at linear order), A should equal A[g + h], where A[-] is just the classical 
functional of the metric which computes the area. Schematically, we have 

(A|h|A) + ---. (3.2) 

h=0 

where A[g] is the classical empty AdS value of the area, equal to (0| A |0) before quan¬ 
tum corrections.^ The result is that the hrst-order corrections to the quantum expec¬ 
tation value for the area can be computed by substituting the expectation value for 
the metric perturbation into the classical area functional, and expanding to first order. 
This is not the case for the higher-order corrections if the metric perturbation behaves 
non-classically, which is as we would expect from a single-particle state with a broad 
wavefunction like |A). 

Furthermore, to hrst order in Gat the expectation value of the metric perturbation 
is obtained easily from the expectation value of the linearized Einstein equations: 

E,4{h)] = (e,„|/!]) = SwGn (t") , (3.3) 

where we have introduced the notation to denote the linearized equation of motion 
for the metric perturbation in the absence of matter (but including the cosmological 
constant).® To summarize, at this order in perturbation theory we can obtain the 
correct backreaction of the matter on the geometry by treating ^ classi¬ 
cal source. From now on we will simplify our notation by writing and 

/^matter \ _ -^matter 

/ih' 

Now it is a simple matter to hnd 6A{T). The hrst-order correction to the extremal 
surface area is obtained by leaving the surface hxed and integrating the linear deviation 
of the area functional over that surface: 

(3.4) 

h=0 

®As pointed out in Ref. [4], quantum corrections to the area term are also important at 0(1). But 
since quantum corrections are already down by a power of N'^, the difference in those corrections 
between the vacuum and first excited states will be further suppressed. 

can still contain contributions from gravitons, and the excited state in question could be a 
pure graviton state. In that case, the metric perturbation would still have to self-consistently satisfy 
the Einstein equations as we have written them to this order in Gat. For example, the expectation 
value of the metric at infinity must encode the total energy of the graviton state. 


6A{T) = e 
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{A\A\A) = {A\A[g + h]\A) = A[g]+ — 
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Bulk Entropy 

The calculation of the bulk entropy is similar to the CFT calculation in Ref. [1]. In 
the low-temperature limit, the reduced density operator for the region S is found by 
taking the trace of (1.5) over the complement of S: 

P^(T) = pI + - pp + ■ ■ • (3.5) 

The second term is just a small perturbation, so we can apply the hrst law of entan¬ 
glement entropy: 

5^buik(T) = -e-^^Tr - p^) log pi] = (3.6) 

where we have denoted by the change in the expectation value of the vacuum 
modular Hamiltonian on E, — logp^. Recall that S is a constant-time surface in the 
AdS-Rindler space constructed in Sec. 2. Much like flat Rindler space, the vacuum state 
of AdS-Rindler is obtained by a 27r tranlsation in Euclidean AdS-Rindler time. So the 
vacuum modular Hamiltonian in AdS-Rindler space is just 27r times time-translation 
Hamiltonian associated to that space, i.e., translation by which is the integral over 
E of one of the components of the bulk stress tensor. Thus we have 

5^buik(T) = 27re-^'^ ^ eT”"**-dE^ (3.7) 

where = dt is the AdS-Rindler Killing vector introduced previously, and dE“ is the 
volume form on E. 


Putting Things Together 

Now we just have to put together (3.4) and (3.7) to get 


^sr°(r) = + 5Sbaii.(r) 

4(jrAr 


(3.8) 

(3.9) 


To see that this is equivalent to in (1.1), we make use of the {d — l)-form 

X constructed out of the metric perturbation h in Ref. [5], which has the (off-shell) 
properties 
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Additionally, the integral of y over the surface B, which from the bulk point of view is 
the part of the boundary of S located at inhnity, can be written as 

[ x = [ d<^-^x6sr[% (3.12) 

Jb J b 

where we still need to identify To do so, recall the holographic dictionary for 

the metric perturbation near the boundary (schematically): 

4 ^( 5 -'')' (3.13) 

where 5 is the vacuum-subtracted expectation value of the stress tensor in the 

CFT. The region B is simple enough that SKb, the vacuum-subtracted expectation 
value of the CFT vacuum modular Hamiltonian on H, can be written as a local integral 

oi5{T^ry 

6Kb = —Tr [{pb — Ps) logp^ = [ (local function of 5 . (3-14) 

Jb 

So, applying (3.13), SKb can be written as a local integral of the asymptotic value of 
the metric perturbation h in the bulk. That local integral is precisely the right-hand 
side of (3.12), which dehnes . Said another way, the integral of y over B is 6Kb 

in the CFT translated into bulk language using the holographic dictionary (3.13). 

If we evalutate (3.12) in the state |A) we will hnd 6Kb in the state |A). By 
including the Boltzmann factor, we obtain the hrst thermal correction to the CFT 
entanglement entropy of H in a way entirely parallel with the bulk calculation in (3.6), 
and in fact this is how (1.1) was obtained by Herzog [1]. In other words, we have 


[ X = = 55^^^(T), (3.15) 

Jb Jb 

where 5S^^'^{T) is the quantity appearing in (1.1). We would like to emphasize that 
this follows directly from the dehnition of and the property (3.12) of y. 

Applying Stokes’ theorem to y yields' 
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The linearized Enstein equations, then give 


iSCFTjj.) ^ ^ I er, 


Jil' ^ flU 
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+ 5^buik(T) = 5^^°'°(T), 
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(3.17) 


^The sign coventions are such that the boundary of S is dX, = B — B. 
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completing the argument. Note that this guarantees that even without 

knowing the explicit form of As an illustration, we go through an example 

in Appendix A where the exact expression in (1.1) is reproduced directly from the bulk. 

4 Discussion 

We have shown that the hrst thermal corrections to CFT entanglement entropy can be 
produced from a purely bulk calculation. There are two main ideas that came together 
to accomplish this. First, as was already noticed in Refs. [5, 7], there is a deep connec¬ 
tion between the linearized Einstein equations in the bulk, the modular Hamiltonian on 
the boundary, and the linearized area functional on the extremal surface. The context 
here is a little different, though. In Refs. [5, 7], the bulk geometry was assumed clas¬ 
sical, and the differential form y was used to derive the linearized Einstein equations 
for that classical metric from the properties of entanglement entropy on the boundary. 
Here we are assuming the Einstein equations, in the form of operator equations inside of 
expectaion values, and using them to learn something about the entanglement entropy. 
The differential form y connects the bulk and the boundary in the same way, but the 
direction of the argument has switched around. 

The second main idea was in our treatment of the Ryu-Takayanagi area functional. 
The metric backreaction was not assumed classical, but still gave a contribution to 
the area of the extremal surface. The idea that the Ryu-Takayanagi term should be 
thought of as an expectation value is very natural, but is difficult to test in a nontrivial 
way. Here we really only used the mild assumption that it behaves as an expectation 
value with respect to thermal averaging and linear quantum corrections. For higher- 
order corrections, it would be important to treat this term in the correct way (which 
likely involves using the generalized entropy in the bulk instead of the area and bulk 
entanglement separately [8]). It would be very interesting if a more nontrivial example 
could be found, such as one where off-diagonal matrix elements of A were important. 
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A Example: Bulk Scalar Field 


In this appendix we will work through an example where the lowest-energy state is that 
of a free scalar field in the bulk, of mass = A{A — d). The lowest energy state is 
the one with zero angular momentum and zero radial quantum number. Let a^, a be 
the creation and annihilation operators for this state, in terms of which the field can 
be written as 

(j) = N cos^ p + ■ ■ • , (A.l) 

where is a normalization constant to be fixed later and the ■ ■ ■ represent cre¬ 
ation/annihilation operators for other modes which do not concern us. In the state 
I A) = |0), the the expectation value of the energy-momentum tensor is spherically 
symmetric and time-independent, and in particular the energy density is 

^matter _ ^2^ cOS^^'^ p. (A.2) 


The coefficient can be computed by demanding that the total energy is equal to 
A. In fact, all we need to know about the energy-momentum tensor is that it is time- 
independent, spherically symmetric, and that the total energy is A. Because of the 
symmetries, the perturbed metric looks like a black hole with a radially-varying mass. 
It is most convenient to add the perturbation to the global metric using the coordinate 
r = tan p, in which case we have to first order 
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It’s straightforward to check that the deviation in the area of B due to this metric 
perturbation is 


5A{T) 


e ^^Vtd -2 
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where we use sin pcos 6 = cos 6q for the surface B. By substituting the expression (2.9) 
for into equation (3.7) for the bulk entropy we find 

^^buik(T) = 27rf2,_2e-^^ / dpddsin-^-^ (A.6) 

is sm do 


9 










We know from the general discnssion of the text that we shonld make nse of Einstein’s 
eqnation to hnd an expression for p in terms of In this case, the rr-component 

of Einstein’s eqnation says 


IottGtv ^ dr 

as is evident from (A.4). Snbstitnting this in for gives 

rr. -/3A / . • rf_2flSinpcos0 - cos^odp 

dh^buikG ) = —^ / dph^sin 9 -—-—. 

sGat Js sint^o dp 

Being clever, we can choose to rewrite the integrand as 
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This allows us to use Green’s theorem: 
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The hrst term 


reproduces the function in (1.1), 

7 -I 


sin p=cos 00 / cos 9 

(A.12) 




while the second term gives us the change in area. 
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(A. 14) 


So we have successfully derived (1.1) from a bulk calculation. 
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